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Optical model potentials for elastic nucleon nucleus scattering are cal- 
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t^j- , culated for a number of target nuclides from a full-folding integral of two 

different realistic target density matrices together with full off-shell nucleon- 
nucleon t-matrices derived from two different Bonn meson exchange mod- 
els. Elastic proton and neutron scattering observables calculated from these 
full-folding optical potentials are compared to those obtained from 'optimum 
factorized' approximations in the energy regime between 65 and 400 MeV 
projectile energy. The optimum factorized form is found to provide a good 
approximation to elastic scattering observables obtained from the full-folding 
optical potentials, although the potentials differ somewhat in the structure of 
their nonlocality. 
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I. INTRODUCTION 



Progress in the rigorous treatment of the multiple scattering of nucleons from nuclei 
has led to the need to study the influence of the full target nuclear density matrix in the 
scattering observables. In first order the spectator expansion of multiple scattering theory 
requires a convolution of the fully off-shell nucleon-nucleon (NN) scattering amplitude with 
the nuclear wave functions of the target. This opens the possibility to assess the influence 
of the target wave functions on elastic proton and neutron scattering observables. 

In its most general form, the first order single scattering optical potential within 
the framework of the spectator expansion is given by the triangle graph shown in Fig. 1. 
Since there is one loop, the evaluation of the graph requires a three-dimensional integration 
involving the fully-off-shell two-nucleon scattering amplitude and a realistic nuclear density 
matrix. Usually, one makes the assumption that the NN amplitude varies slowly as a 
function of its arguments compared to the nuclear density matrix. This corresponds to 
the argument that the range of the NN force is small compared to the size of the nucleus 
and leads to the approximate nonrelativistic form t(q)p(q) for the first-order nucleon- nucleus 
optical potential. Full-folding calculations, avoiding this approximation, have already been 
performed by several groups based on either the KMT approach 0|| or a g-matrix approach 
using various models for the off-shell density matrix as well as different models for the 
NN amplitudes. In general, this work indicates that an improved treatment of the off-shell 
structure of the optical potential improves the description of the observables. 

Our approach to elastic scattering from nuclei is based on the spectator expansion 
of multiple scattering theory Here the first order term involves two-body interactions 

between the projectile and one of the target nucleons. Due to the many-body nature of 
the free propagator for the projectile+target nucleon system, a theoretical treatment of this 
many-body propagator as affected by the residual target nucleus is included. The calculation 
of the optical potential presented in this paper relies on two basic inputs. One is the fully off- 
shell NN t-matrix, which represents the current understanding of the nuclear force, and the 
other is the nuclear wave function of the target, representing the best understanding of the 
ground state of the target nucleus. To account for the modifications of the free propagator 
inside the nucleus, the same mean field potentials are used from which the ground state 
wave functions are derived. There are no adjustable parameters present in this calculation. 

The motivation for ongoing work on this topic is twofold. First, elastic and in- 
elastic nucleon-nucleus scattering provide an important and sensitive test for theoretical 
corrections at the first-order level of the optical potential (e.g. as given by possibly gen- 
uine modifications of the NN interaction in the nuclear environment and off-shell effects). 
Rigorous microscopic calculations are required for discerning these effects. Second, a bet- 
ter understanding of the theoretical details of the optical potential are needed to construct 
realistic and physically sound wave functions representing continuum nucleons in the in- 
terior of the nucleus. These wave functions will become vital for future theoretical needs 
in high-energy coincidence experiments (at TJNAF e.g.), inelastic scattering studies, and 
for understanding the reactions in heavy ion experiments involving nuclei far from the drip 
lines. Here experiments involving the scattering of exotic nuclei from single nucleon targets 
should should benefit from full-folding type calculations in order to test the predicted density 
distributions of halo-like nuclei. 
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The structure of the paper is as follows. In Section II we review the relevant 
formalism for the single-scattering optical potential and introduce the full-folding procedure 
as used in our calculations. In Section III we discuss the model densities employed and 
describe the calculations of the full folding optical potentials. Elastic scattering results for 
neutron and proton scattering form a variety of nuclei in the energy regime between 65 and 
400 MeV are discussed in Section IV. Our conclusions are presented in Section V. 



II. THE FULL-FOLDING OPTICAL POTENTIAL 



The standard approach to elastic scattering of a strongly interacting projectile 
from a target of A particles is the separation of the Lippmann-Schwinger equation for the 
transition amplitude 

T — V + VG (E)T (2.1) 
into two parts, namely an integral equation for T: 

T = U + UG (E)PT, (2.2) 
where U is the optical potential operator and defined by a second integral equation 

U = V + VG (E)QU. (2.3) 

In the above equations the operator V represents the external interactions between the 
projectile and the target nucleons. Therefore the Hamiltonian for the (A + 1) particle 
system is given by 

H = H + V. (2.4) 

The free propagator G (E) for the projectile-target system is given by 

G (E) = (E-H + ie)-\ (2.5) 

The potential operator V = Yli=\ v m consists of the two-body potential v 0i acting between 
the projectile and the zth target nucleon. The operators P and Q are projection operators, 
P + Q = 1, and P is defined such that Eq. (|2.2|) is solvable. In this case, P is conventionally 
taken to project onto the elastic channel, such that [Go? P] — 0. The free Hamiltonian is 
given by 

H = h + H A (2.6) 

where h is the kinetic energy operator for the projectile and Ha stands for the target 
Hamiltonian. Thus the projector P can be defined as 



p _ |$a)($a| 

($a|$a) 

where |<E>a) corresponds to the ground state of the target and fulfills 



(2.7) 



3 



H A \$ A ) = E A \® A ) (2.8) 

With these definitions the transition operator for elastic scattering can be defined as 
T e i = PTP, in which case Eqs. fl2.2p becomes 



T el = PUP + PUPG {E)T el . (2.9) 

The fundamental idea of the spectator expansion for the optical potential is an 
ordering of the scattering process according to the number of active target nucleons interact- 
ing directly with the projectile. The first order term involves two-body interactions between 
the projectile and one of the target nucleons, i.e. 

t/ = X>, (2.10) 

i=l 

where the operator T; is derived to be 

n = v 0i + v 0i Go(E)QTi 

= v 0i + VQiGoffln - v 0i G (E)Pr, l (2.11) 
= ti- f l G (E)Pr, l . 

For elastic scattering only PTjP, or equivalently (^aI^I^a) need to be considered, 

($aM$a) = (Q A \Ti\$A) - (QA\Ti\$ A )TB tp I u I ■ ( $ aN$a), (2.12) 

(E - E A ) —ho+ie 

where fj is defined as the solution of 

n = v 0i + v 0i G (E)fi. (2.13) 

It should be noted that Eqs. ( ^.3[ ) to ( |2.13| ) all follow in a straightforward derivation and 
correspond to the first order Watson scattering expansion 0. If the projectile — target 
nucleon interaction is assumed to be the same for all target nucleons and if isospin effects 
are neglected then the KMT scattering integral equation [[Kj can be directly derived from 
the first order Watson scattering expansion. 

Since Eq. ( P-12j ) is a simple one-body integral equation, the principal problem 
is to find a solution of Eq. (|2.13|) , which has a many-body character due to Gq(E) = 
(E — h — H A + ie)' 1 . If the propagator Gq(E) is expanded in the the spirit of the spectator 
expansion within a single particle description, one obtains in first order [fUJIl}] 

Gi{E) = [{E - E') -h -hi-Wi + ie}-\ (2.14) 

where hi is the kinetic energy of the ith. target particle and Wi = J2j^i v ij- The quantity Wj 
represents the force acting between the struck target nucleon and the residual (A-l) nucleus. 
Then the operator fj of Eq. (|2.13| ) can be written as 

n = v 0i + voiGi(E)fi 

= t Q i + t oi g i W i G i (E)f i . (2.15) 
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Here the operators toi and are defined as 

tot = v oi + v 0i git oi (2.16) 

and 

g i = [(E-E i )-h Q -h i + ie]- 1 . (2.17) 
The operator t 0i can be identified with the free NN t-matrix, and in lowest order the operator 



fi of Eq. (|2.15| ) is given by fj ~ toi- From now on we consider for clarity in presentation only 
this case. 

The matrix element (^aI^I^a) given in Eq. ( |2.12j ) represents the full- folding op- 
tical potential and is given explicitly as 

(k'|C/|k) = (k'^l J2 r a \k*A), (2.18) 

a=p,n 

where a represents the sum over the target protons and neutrons. Since (k'|{7|k) is the 
solution of the sum of the one-body integral equations represented by Eq. ( |2.12| ), it is 
sufficient to consider the driving term 

(k'\U\k) = (k'$4| ]T ^|k$4>, (2.19) 

a=p,n 

where f a is given by Eq. ( |2.13| ). 

Inserting a complete set of momenta for the struck target nucleon before and after 
the collision Eq. fl2.19p reads 

U (k, k') = £ / d 3 p'd 3 p (k'p' | f a (e) | kp) p a (p' + p p + ^ 5 3 (k' + p' - k - p). (2.20) 

The momenta k' and k are the final and initial momenta of the projectile in the frame of zero 
total nucleon- nucleus momentum. The structure of Eq. ( [2.20 ) is represented graphically by 



Fig 1, which also illustrates the momenta p' and p. The proton and neutron density matrices 
are given by p a . By evaluating the ^-function and introducing the variables q = k' — k, 
K = |(k + k') and p = |(p' + p) we obtain 



£>(q,K)= £ J d 3 p ^k',p-iq| f Q (e) |k,p+iq 



a=p,n 

( A -la K A-la K\ 
"H P "~2 + ^ P+ ~2 + ^)' (2 ' 21) 

A change of the integration variable from p to P = p + accounting for the recoil of the 
nucleus, gives [O 



f>(q,K)= £ /d 3 P (k',P-f-f |r a ( e )|k,P + |-^ 



a=p,n 



A - 1 q _ A - 1 q\ 



The NN amplitude f a in Eq. (|2.22|) is evaluated in the zero momentum frame of the nucleon- 
nucleus system. The relationship to the corresponding matrix element evaluated in the zero 
momentum frame of the two nucleons is given by 




rKP,q,K)<r,-/C'| f Q (e) \Kl,-K) NN , (2.23) 



NA 



where KJ = W - (P 



q 

2 



)) and K 



(k-(P 



)) are the nonrelativistic final 



and initial nuclear momenta in the zero momentum frame of the NN system. The factor 
T) (P, q, K) is the M0ller factor for the frame transformation [14]] and is given by 



7?(P,q,K) 



E N (JC')E N (-IC')E N (IC)E N (-IC) 



E N (k')E N (P-2-K)E N (k)E 



;p + f 



K 

A > 



[2.24) 



where E^(k) is the relativistic kinetic energy of a nucleon of momentum k. This factor 
imposes the Lorentz invariance of the flux. With this frame transformation taken into 
account, the full-folding optical potential of Eq. ( |2.22|) can be written as 



U(q,K)= £ [d 3 Pr,(P,q,K)f a U 

a=p,n 



A 



K-P);e 



A - 1 q „ A - la 
Pa P — x,P + 



A 2 



A 2 



(2.25) 



Here the arguments of f a are q = k' - k = JC' - K and \{KJ + K.) = |(^K - P). 

The two-nucleon amplitude f a is calculated from the free NN t-matrix according 
to Eqs. ( |2.15| ) and ( 2.16| ) at an appropriate energy e. In principle, this energy should be the 
beam energy minus the kinetic energy of the center of mass of the interacting pair less the 
binding energy of the struck particle. Following this argument, e should be coupled to the 
integration variable P. The full- folding calculations of Refs. |4|,[| are carried out along this 
vain. For our calculations we take a different approach, namely we fix e at the two-body 
center-of-mass (cm.) energy corresponding to free NN scattering at the beam energy so that 
the same laboratory energy applies to the two-body and nuclear scattering. This approach 
has been applied in earlier work || and is also used in the work of the Surrey Group |J. 



III. MODELS FOR THE OFF-SHELL DENSITY 



The evaluation of the full-folding optical potential as given in Eq. Q2.25 ) requires 
a nuclear density matrix, which in a single particle picture is given as 

P«(P / ,P)=E*U^)*«,*(P) 

i 

Here \I/ Qj j(p) are the wave functions describing the single particle nuclear ground state. The 
index a stands for protons and neutrons, respectively, and the total nuclear ground state is 
given by the sum of the two. In order to achieve consistency with our formulation of incorpo- 
rating effects of the 'nuclear medium' on the scattering process we choose as model density 
matrices the ones from which the nuclear mean fields Wi are derived (see e.g. Eq. (|2.15|) ). 
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The models used are a non-relativistic reduction of a Dirac-Hartree calculation ]TJJ and a 



non-relativistic Hartree-Fock-Bogolyubov (HFB) structure calculation [ 17,18] . The Dirac- 
Hartree calculation is a spherical solution of the one-body Dirac equation assuming a scalar 
potential and the time component of a vector potential field. The nonrelativistic HFB mi- 
croscopic nuclear structure calculation uses the parameterized effective finite-range, density 
dependent Gogny D1S effective NN interaction. The parameter of the Gogny D1S interac- 
tion are fitted to a certain set of stable nuclei. For this case an axial harmonic oscillator 
basis is used. 

The details of the Dirac-Hartree (DH) calculation leading to the density matrices 
employed in our calculations are given below. The wave function \&i(r) is a solution of the 
one-body Dirac-Hartree equation and given by [|19| 



n,v,m,t z 



j- G*s,n,i/(r) j 

•jWMi 



(3.2) 



Here t z stands for the z component of the isospin and n for the principal quantum number. 



The phase convention is taken from Ref. |T{|. During this derivation we prefer to omit the 
index a. The spherical harmonics are determined by <fi vm which is defined as 



1 1 

w= zl ( im i~ m s\hj m ) Y r i (tyxr. 



(3.3) 



where Y^iQ) is a spherical harmonic and Xm a a Pauli spinor. The quantum number v 
uniquely defines j and I as 



J 



1 
2 



I 



is, v > 
■(y + 1), z/<0 



(3.4) 



We used the code Timora JT9|] and the parameter sets given therein to generate the functions 
G and F given in Eq.(|3.2|) for the nuclei studied in this paper. 

Under the assumption of orthogonal single particle states the density matrix is given in 
coordinate space by 



p(r',r 
G 



numt. 



E 



/ : ,n,ui r ) Gt z n ,u( r ) sr^ i i i\ , / \ . F z ,n,u( r ) F z ,n,u( r ) sr^. < , n< , \ 

2^ <Pvm{r )4>vm{r) H 2^ (P-vmKr )(p~vm{r) 



. (3.5) 



Here we should point out that in order to obtain a density matrix which we can apply 
in our formulation of the optical potential, we have a 1-operator between the Dirac wave 
functions ^ num t z , and then treat p(r', r) as nonrelativistic density matrix. The orthogonality 
of the spin states leads to 5 ms , m ' s and thus to mi = ray. Taking advantage of the symmetry 
properties of the Clebsch-Gordon coefficients leads to 



P( r » = E 



G tz ,n, v (r') Gt z ,n,v(r) F tz ^ v (r>) F u ^ v {r) 



2j_+l 
21 



J2Yr(r')Yr( 



(3.6) 



The calculation of the full-folding optical potential J7(q, K) requires the nuclear 
density matrix in momentum space. Thus we need to double Fourier transform p(r', r) to 
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obtain the density p(p', p) in the rest frame of the nucleus. This frame is characterized by 
the momenta p and p' and the density matrix is obtained by 



Pc 



P) = g^3 / d 3 r'e- ir '-v' J d 3 re^p a (r', r), (3.7) 



where we again indicate with the index a that we have to obtain the density matrix for 
protons as well as neutrons. Using the standard expansion of a plane wave, the angular 
integration in Eq. (|3.7|) can be easily carried out, and we obtain for the density matrix 



P Q (P',P) = ^E( 2J+1 )E P KC0S^ 



™ J 



j "*V j^WO Jdrr Jl{ pr)F a ^ A r) + 

dr'r' ji(p'r')G^t z ,nA r ') / dr r 3i(P r )G a ,t z ,nA r ) ]■ ( 3 - 8 ) 



The density matrix p a (p',p) given in Eqs. (|3.7|) or (3.8) is defined in the rest 
frame of the nucleus. In order to apply p a {p', p) in our calculation of the full- folding optical 
potential for nucleon-nucleus scattering, we have to evaluate the function at the correspond- 



ing momenta in the nucleon-nucleus frame. This is facilitated by variable transformations 
J and p' = p' - 



p = p — x and p' = p' — 7, which takes into account recoil. As an aside, not including 



A' 



recoil would mean the transformation p' = p' 

For the calculation of the density matrices derived from a non-relativistic Hartree- 
Fock Bogolyubov (HFB) calculation based on the Gogny-DlS NN interaction we employ 
essentially the same procedure as described above. The wave functions are created in r 



space by a code provided by Berger |l7j and are represented in a axially symmetric harmonic 
oscillator basis. A double Fourier transform is then performed using the oscillator basis and 
summing over all harmonic oscillator quantum numbers. This choice of basis takes advantage 
of the fact that the Fourier transform of a harmonic oscillator is again a harmonic oscillator. 
The density matrix is given by 



p4p >pj 



£P M V1'(P>*(P), (3.9) 



where the indices i, %' count the harmonic oscillator basis states and p 1 ' 1 ' is the density matrix 
in the oscillator basis. Again, the index a distinguishes between protons and neutrons. The 
basis states are explicitly given by 

^(p) = £ Am,(/3,7) e"^ 2 H t ((3,p z ) e~^ 2 (l,Pr)e vmd . (3.10) 

m 

Here p z is the projection of the momentum along the z-axis and p r the radial momentum, 
P, 7 are harmonic oscillator constants. Hi{(3^p z ) are the Hermite polynomials and L| m '(7,p r ) 
the Laguerre polynomials. The size of the harmonic oscillator basis used depends on the size 
of the nucleus, e.g. the size of the basis for 16 is 12 shells whereas for 90 Zr it is 16 shells. 
It should be noted that the indices i and %' are not independent. The size of the basis sets 
needed makes the calculation of p a (p',p) quite lengthy, especially for heavier nuclei. 
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In order to calculate the optical potential U (q, K) as given in Eq. (|2.25|) , we need 



the density matrix as function of the momentum transfer q and P = p + as indicated in 
Eq. (|2.22| ). In these variables the density matrix is related to the density profile p a (<?) of the 
nucleus by 

p a (g) = /rf 3 Pp a (p-^f,P + ^f). (3.11) 

The normalization is chosen such that p a (q = 0) = Z or N, the number of protons or 
neutrons, respectively. 

In practice we used the relation given in Eq. ( |3.11|) for testing our numerical in- 
tegration schemes with the simple harmonic oscillator density given in Ref. ||. In order 
to determine how well the two model density matrices presented here describe the experi- 
mentally determined proton distribution, we calculate the proton density profiles p p (q) for 
both the DH and the HFB models for each nucleus we consider. In the following we want to 
discuss two cases, namely 16 and 90 Zr. In Fig. 2 we compare the density profiles calculated 
from the DH and HFB models to the experimental proton distribution PU| . Overall the DH 
profile follows the experimental distribution closer than the HFB profile. The HFB profile is 
shifted to larger momenta indicating that the HFB model slightly underpredicts the radius 
of the proton distribution of 16 O. This feature will be visible in the proton-scattering observ- 
ables for 16 O calculated with the HFB model. In the close-up of the minimum of the density 
profile it can be seen that both model densities slightly deviate from the experimental profile. 
In Fig. 3 we carry out the same comparison for a heavier nucleus, 90 Zr. Here both model 
densities follow the experimental proton distribution |2U] very closely. The close-up of the 



minimum reveals that the HFB profile deviates only slightly from the experimental profile. 
This is a general trend, the heavier nuclei are better described by the model profiles. In fact, 
the proton distribution of 16 represents the worst case of disagreement of the model profiles 
with the experimental profiles. This is understandable since the HFB model is known to 
provide a better representation of the larger nuclei. 



IV. RESULTS AND DISCUSSION 

A. Details of the Calculation 

In this paper the study of the elastic scattering of neutrons and protons from spin- 
zero target nuclei at energies that range from 65 to 400 MeV (incident projectile energy) is 
strictly first order in the spectator expansion. Here the connection to the propagator Gq{E) 
due to the coupling of the initially struck target nucleon to the residual target is considered 
to be first-order. The full-folding optical potential is calculated as outlined in Section II, 
specifically as given in Eq. ( |2.25| ), using the model densities described in Section III. The 
calculations for scattering at energies smaller than 200 MeV take into account the coupling 
of the struck target nucleon to the residual nucleus via the mean field potential Wi, which 
is chosen to be consistent with the model density employed. Details of this procedure are 
given in Refs. P|JTT[]. Calculations using the Dirac-Hartree densities (and the corresponding 
potential Wi) are labeled DH, while those using the Hartree-Fock-Bogolyubov densities (and 
corresponding potentials Wi) are labeled HFB. 
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The convolution of the fully off-shell density matrix p a with the fully off-shell NN t- 
matrix, and the M0ller frame transformation factor i](P, q, K) as given in Eq. ( p. 24 ) is carried 



out in three dimensions without partial wave decomposition and the integration is performed 
using Monte Carlo integration techniques. Our algorithm uses Quasi-Random numbers 



2jJ , together with importance sampling, which according to our tests has the advantage of 
needing significantly fewer integration points than algorithms based on conventional 'random 
number' generators or Gauss-Legendre integration to obtain the same accuracy. Quasi 
random numbers provide a 'uniform' random distribution over the integration space. 

Aside from the density matrices, the fully off-shell NN t-matrix is another cru- 
cial ingredient in the calculation of C/(q, K). The calculations presented use the free NN 
interaction based upon the full Bonn potential ||22|| . This interaction includes the effects of 



relativistic kinematics, retarded meson propagators as given by time-ordered perturbation 
theory, and iterative and crossed meson-exchanges with NN, NA, and AA intermediate 
states. For the calculations involving projectile energies greater then 300 MeV we employ 



an extension of the Bonn model above pion-production threshold JT^|. In this model pion 
production is described through the decay of the delta isobar with a width obtained con- 
sistently from the imaginary part of the one-pion loop diagram for the delta self-energy. It 
is also to be understood that we perform all spin summations in obtaining £/(q, K). This 
reduces the required NN t-matrix elements to a spin independent component (correspond- 
ing to the Wolfenstein amplitude A) and a spin-orbit component (corresponding to the 
Wolfenstein amplitude C). Since we are assuming that we have spin saturated nuclei, the 
components of the NN t-matrix depending on the spin of the struck nucleon vanish. For the 
proton nucleus scattering calculations the Coulomb interaction between the projectile and 



the target is included using the exact formulation described in Ref. |27] . 



A common approximation to the full-folding expression of Eq. ( |2.25| ), which still 
preserves the non-local character of the NN t-matrix, is obtained as follows. If one observes 
that the nuclear size is significantly larger than the range of the NN interaction, the am- 
plitude f a is expected to be the slower varying quantity in the integral of Eq. ( |2.25|) . This 



argues for the method of optimum factorization [0,0] which proceeds via an expansion of 
f a (including the factor ?](q, K, P)) in P about a fixed value Po- The reference momentum 
Po is determined by requiring that the contribution of the first derivative term be mini- 
mized. In the elastic scattering case this contribution can be made to vanish if Po is chosen 
to be zero. For further details we refer to Ref After the integration over the density 
matrix to produce the diagonal density profile p a (q) (Eq. ( pj.ll| ) the 'optimum factorized' or 
'off-shell tp' form of the optical potential is given by 

U fac (q, K) = £ r)(q, K) f a (q, 4t^K, e) p Q (q) (4.1) 

a=p,n \ ^^i ' 

Here the non-local character of the optical potential is solely determined by the off-shell NN 
t-matrix. For harmonic oscillator model densities it has been shown for light nuclei that 
the optimum factorized form represents the non-local character of U (q, K) qualitatively 
when applied within the KMT formalism to first order at intermediate energies. When 
comparing elastic scattering observables obtained from full-folding optical potentials to those 
obtained from 'off-shell tp' optical potentials, the scope is two-fold. First, we employ here 
realistic models of the nuclear density for light as well as heavy nuclei. Second we extend this 
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comparison toward energies below 100 MeV where it could be expected that the nucleon- 
nucleus scattering calculation samples the optical potential further off-shell and thus the 
optimum factorized form may not be as good an approximation. 



B. Elastic Scattering Results 

Elastic scattering calculations from several spherical nuclei are carried out at a 
variety of energies between 65 and 400 MeV to allow for comparisons between results ob- 
tained from the full-folding optical potentials with those arising from the factorized off-shell 
'tp' approximation. 

The scattering observables for elastic proton scattering from 16 O are displayed 
in Fig. 4. The solid line represents the calculation with the full-folding optical potential 
based on the DH density and the Bonn t-matrix defined above pion-production threshold, 
while the dashed line represents the optimum factorized form as defined in Eq. Q4.1| ). Both 
calculations are based on the free NN t-matrix. Since the two calculations give very similar 
results, it can be concluded that the bulk of the non-locality of the optical potential, which 
affects the elastic scattering observables, must come from the off-shell structure of the NN 
t-matrix. The off-shell structure of the nuclear density matrix plays an insignificant role for 
elastic scattering observables at these high energies. A similar conclusion was already drawn 
in Ref. || and is here confirmed using realistic densities. In order to illustrate the effect 
of the different density profiles for the DH and the HFB models (as shown in Fig. 2) on 
the elastic observables, we display two calculations based on the factorized optical potential 
for the DH (solid line) and the HFB model (dashed line) in Fig. 5. As already discussed 
in Section III, especially in the case of 16 0, the HFB density profile is shifted to larger 
momenta compared to the DH profile. This translates directly into a slight shift of the 
first minimum of the differential cross-section to larger angles and a slightly smaller angular 
spacing of the diffraction minima. We carried out similar comparisons of full-folding and 
optimum factorized optical potentials for heavier nuclei, but there the disagreement between 
the density profiles of the two models is much smaller then for 16 O and consequently the 
prediction of the observables are very similar. 

Another effect worthwhile to study in this context is the influence of the M0ller 



factor [II], which takes into account the transformation of the NN t-matrix evaluated in the 
NN cm. frame to the zero momentum frame of the nucleon-nucleus system. This frame 
transformation can be viewed as a relativistic effect and its importance should increase with 
higher scattering energies. For these reasons we want to consider it's effect on the elastic 
scattering observables for proton scattering from 16 at 500 MeV (Fig. 6). The M0ller factor 



r](P, q, K) as given in Eq. ( |2.24|) is a function of three vector momenta and is part of the full- 



folding integral of Eq. ( 2.25|) . The solid line in Fig. 6 represents the calculation of £/(q, K) 



as given in Eq. ( |2.25| ). In the spirit of the optimum factorized approximation ?7(P,q, K) 
can be expanded around a fixed value P (here P o =0), thus becoming independent of the 
integration variable P. This expansion corresponds to considering r/(q, K) at a fixed angle 
between q and K, specifically here = 90°. The dashed line therefore in Fig. 6 corresponds 
to evaluating the optical potential according to 
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t/(q,K)= •£ ),(q.K)e.90» /i s P f„ (q, - 



a=p,n 



A 



( A - 1 q A - 1 q\ , . 

The dashed and solid lines in Fig. 6 are almost indistinguishable. This infers the conclusion 
that ?](q, K) e =9o° is a very good representation of the exact expression given in Eq. fl2.24| ). 



In order to illustrate the total effect due to the inclusion of the M0ller factor, the dotted 
line in Fig. 6 represents a calculation with 77 (q, K) set to one in Eq. ( |4.2| ). 

At energies below 200 MeV, calculations of elastic observables not only incorporate 
the effects of the nuclear structure models within the full-folding procedure but also via the 
mean field force ( given by the structure model) which couples the struck target nucleon 
to the residual nucleus. Thus it is hoped that the influence of different structure models 
on the elastic observables is observable. In Fig. 7 we display the elastic observables for 
proton scattering from 40 Ca at 100 MeV laboratory energy employing the DH model for the 
density as well as the mean field force W{. In Fig. 8 the corresponding calculation is done 
using the HFB model. In both figures the solid line represents the full-folding calculation, 
and the dashed line the factorized off-shell 'tp' approximation. All calculations contain the 
modification due to the mean field W%. The off-shell structure of the nuclear density matrix 
in the full-folding procedure has at this lower energy a slightly larger effect on the spin 
observables then at higher energies. In addition, the angular distribution of the differential 
cross section diffracts at slightly larger angles in the full-folding calculations compared to 
those based on the factorized form. This trend is also be observed for the heavier nuclei 90 Zr 
and 208 Pb (Figs. 9 and 10). 

The elastic observables for proton scattering from 90 Zr at 80 MeV are displayed 
in Fig. 9. Here the difference between the two model densities employed is almost negligible 
for ^ and A y . Only for the spin rotation function the difference given by using two different 
structure models is at higher angles as large as the effect of using the factorized approxima- 
tion. This result, namely that the observables predicted by the two different models under 
consideration are so similar, is not surprising in that both models predict an almost identical 
density profile (Fig. 3). The effect of the off-shell structure of the nuclear density matrix 
is relatively small as the comparison between the full-folding (solid) line and corresponding 
factorized (dashed) calculation shows. In the case of proton scattering from 208 Pb at 65 
MeV (Fig. 10) a comparison between a full-folding calculation and the factorized approxi- 
mation reveals the same trends as the observables for the lighter nuclei. At large angles the 
full-folding calculation falls below the one given by the factorized form, and in this case also 
below the data. The inclusion of the off-shell structure of the nuclear density matrix makes 
the nucleus appear slightly larger, which becomes apparent in the shifted diffraction pattern 
of — 

It is difficult to extract properties of nonlocal potentials from elastic scattering 
observables. Nonlocal effects are presumably more important in inelastic processes which 
depend on the nucleon-nucleus interaction such as quasielastic electron scattering reactions. 
In order to gain more insight into the difference between a full-folding optical potential and 
the factorized off-shell 'tp' approximation to this potential, we plot in Figs. 11 ( 40 Ca) and 
14 ( 208 Pb) the real and imaginary parts of the on-shell value of t/(q, K) as a function of the 
orbital angular momentum L. We separate the cases J = L + ^ and J = L — | to isolate the 
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effect of the spin-orbit force. As is seen in both figures, the full-folding (solid lines) and the 
factorized (dashed lines) on-shell values of the imaginary parts of the optical potential are 
quite similar. In both cases the real part of the on-shell values of the optical potentials exhibit 
an increasing suppression for smaller L as the nonlocal effects of the density matrix as well 
as the NN t-matrix are treated more adequately in the full-folding procedure. It should be 
emphasized again that Figs. 11 and 14 only show the value of C/(q, K) fulfilling the on-shell 
condition q • K = and q 2 + 4K 2 = 4ko 2 with ko being the on-shell relative momentum 
for proton-nucleus (NA) scattering and do not display any off-shell behavior inherent in 
the potentials. After iteration to obtain the Watson optical potential (Eq. |2.11| ) and then 
in the integral equation (Eq. |2.9| ), the on-shell elements of the NA t-matrix display much 
smaller differences between the full-folding and it's factorized, on-shell '£p' approximation. 
For 208 Pb the differences in the real parts are nearly insignificant (Fig. 15), whereas for 
40 Ca a slight suppression of the real part for small L remains for the full-folding calculation 
compared to the factorized approximation. However, the differences occur mainly for smaller 
L where the imaginary part of the potentials as well as the t-matrices are relatively large. 
Because the absorption is significant for these low partial waves, the elastic observables are 
not particularly sensitive to the differences in the real parts as displayed in Figs. 13 and 16. 

In Fig. 17 total neutron cross section data for 12 C, 16 0, 28 Si, 40 Ca, 90 Zr, and 
208 Pb are shown along with various calculations of a tot (E) at a number of energies. Because 
the data are so extensive, the 'usual' procedure has been reversed in the plotting of these 
cross sections so that the data are represented by dotted curves, and the discrete points 
correspond to calculated results. The solid diamonds represent the full-folding calculations 
as described in Section II. All calculations are based on a DH model for the nuclear density. 
For energies < 200 MeV the modification of the free propagator through the DH mean field 



is included as described in Ref. [IT|. It has been shown in Ref. ][]] that for higher energies 
this modification of the free propagator becomes negligible. The open circles represent 
calculations based on the factorized, off-shell 'tp' form using the same NN t-matrix. A 
general trend to be observed in Fig. 17 is the slightly lower value of <Jtot{E) obtained from 
a full-folding calculation compared to the factorized approximation. This trend is almost 
independent of the energy and the nucleus under consideration and is consistent with the 
observation that full-folding calculations of the differential cross sections fell slightly below 
the values given by a factorized calculation. 

At this point it is worthwhile to investigate whether the interactions of the pro- 
jectile with the target nucleus is uniformly distributed over the entire nucleus or if specific 
regions of the nucleus play a more dominant role in the scattering process at intermediate 
energies. For our study we chose neutron scattering at 200 MeV projectile energy and con- 
sider contributions from specific shells to the total cross sections. We employ the DH model 
for the density and remove outer shells of protons as well as neutrons. Then we recalculate 
the scattering from the remaining 'inner core', which is chosen to be doubly magic, so that 
it is bound. The results of this procedure for 16 0, 40 Ca, and 208 Pb are given in Table I. As 
a technical detail, when calculating the scattering in these tests we treated the targets as 
being infinitely heavy to exclude recoil effects, which would be larger for smaller cores. In 
order to give an estimate of the size of the recoil effect on the total cross section we give the 
values of Otot calculated with and without recoil in Table I. The values for the total cross 
section for neutron scattering for 'inner cores' of 100, 40, 16, and 4 nucleons are given as 
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entries of the corresponding nuclei. The entries in Table I marked 'n.b.' indicate that e.g. 
in the case of 208 Pb the DH calculation with only 8 neutrons and 8 protons did not result in 
a bound system using the parameters of 208 Pb given for those 16 nuclei. The calculated rms 
radii for the 'inner cores' under consideration for 16 0, 40 Ca, and 208 Pb are listed in Table II. 
This table also contains the rms radii for the proton and neutron distributions for the above 
mentioned nuclei as given by the DH model. Columns 3 and 4 of Table II compare the 
percentage of the volume filled by the 'inner core' if either the corresponding rms radius is 
used (column 3) or the radius is taken to be proportional to A 3 (column 4). The percentage 
Of the calculated total cross section contribution from the inner core nucleus is given in 
column 5. The numbers suggest that the nucleons in the interior of the nucleus contribute 
to the total cross sections with a percentage slightly larger then the volume they occupy 
when the volume is based on the crude estimate r ~ As. This leads to the conclusion that all 
nucleons in the nucleus almost equally contribute to the scattering process. We performed 
a similar study at 100 MeV and 500 MeV projectile energy and did not find any significant 
deviations from the ratios a tot (core) / a tot as given in Table II at 200 MeV. 

V. CONCLUSION 

We have calculated the full-folding integral for the first-order optical potential 
within the framework of the spectator expansion of multiple scattering theory. These optical 
potentials are based on realistic models for the nuclear density matrix, namely a Dirac- 
Hartree and a Hartree-Fock-Bogolyubov model along with the full Bonn meson exchange 
model for the NN t-matrix. Recoil and frame transformation factors are implemented in 
the calculation in their complete form. We calculated elastic scattering observables for a 
variety of light and heavy nuclei at projectile energies from 65 to ~400 MeV laboratory 
energy. At energies below 200 MeV we included the modification of the free propagator 
due to the coupling of the struck target nucleon to the residual nucleus via the same mean 
field used to model the effect of the nuclear medium is intrinsically consistent with the 
nuclear structure. The predictions from these rigorous calculations of elastic nucleon nucleus 
observables provide excellent agreement with the experimental data in this energy regime. 

We tested the validity of the factorized off-shell 'tp' approximation in the energy 
regime between 65 and 400 MeV and found that this approximation, which only retains the 
non-locality given through the NN t-matrix, is even at lower energies a very good repre- 
sentation of the full-folding calculation as far as the elastic nucleon-nucleus observables are 
concerned. Differences between the factorized approximation and the full calculation of the 
optical potential are present predominantly in lower partial waves. However due to the cu- 
mulative effect of many partial waves the elastic observables do not reflect these differences. 
It should be noted that in e.g. inelastic scattering of nucleons from nuclei or quasielastic 
electron scattering those differences between full-folding calculations and the corresponding 
factorized approximation may become more significant. We also studied the contribution of 
the interior structure of the nucleus to the total cross section and find that all nucleons in 
the nucleus contribute almost uniformally to the scattering process. 
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TABLES 



TABLE I. Total cross sections for neutron scattering from 6 0, Ca, and 08 Pb as well as from 
inner shells of those nuclei. The entries printed in boldface are the ones for which the rms-radii 
and ratios are calculated in Table II. 





vtot[b] 


crtot[b] 


&tot[b] 


crtot [b] 


&tot[b] 


&tot[b] 






no recoil 


core of 4 


core of 16 


core of 40 


core of 100 


16Q 


.423 


.419 


.120 


.419 






40 Ca 


.925 


.921 


n.b. 


.419 


.921 




208p b 


3.38 


3.37 


n.b. 


n.b. 


.960 


2.01 



TABLE II. Rms-radii for the proton and neutron distributions of 16 0, 40 Ca, and 208 Pb as well 
as from inner shells (cores). The last three columns give the ratios of the volumes of the cores to 
the total nucleus as well as the ratios of the calculated total neutron cross sections. The numbers 
used to determine the latter are the ones printed in bold in Table I. 





rms-radius(full) [fm] 


rms-radius(core) [fm] 






atot(core) 
crtot 




(proton, neutron) 


[core] : (proton,neutron) 


{vTflS cor e^) 

(rms 3 ) 


A 


16Q 


(2.63, 2.60) 


[4]:(1.96, 1.95) 


42% 


25% 


29% 


40 Ca 


(3.39, 3.33) 


[16]:(2.63, 2.60) 


48% 


40% 


45% 


208p b 


(5.40, 5.67) 


[40]:(3.72, 4.91) 


47% 


19% 


29% 
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FIGURES 



FIG. 1. Diagram for the optical potential matrix element for the single-scattering term. 



FIG. 2. Comparison of the experimental proton density profile p p (q) for 16 |2(| (solid line) 
with the calculated proton density profiles from the DH model (dash-dotted line) and the HFB 
model (dashed line). 



FIG. 3. Comparison of the experimental proton density profile p p (q) for 90 Zr ]^(J (solid line) 
with the calculated proton density profiles from the DH model (dash-dotted line) and the HFB 
model (dashed line). 

FIG. 4. The angular distribution of the differential cross-section (^), analyzing power (A y ) and 
spin rotation function (Q) are shown for elastic proton scattering from 16 O at 400 MeV laboratory 
energy. The solid line represents the calculation performed with a first-order full-folding optical 
potential based on the DH density [16| and the Bonn model D52 [jlf|. The dashed line represent 



the calculation using the factorized, off-shell 'tp' approximation to this optical potential. The data 
are taken from Ref. [SO]. 



FIG. 5. Same as Fig. 4, except that all calculations are based on the factorized, off-shell 'ip' 



approximation. The solid line represents the calculation using the DH [16| density profile, whereas 
the dashed line uses the HFB 17] density profile. 



FIG. 6. The angular distribution of the differential cross-section (^), analyzing power (A y ) 
and spin rotation function (Q) are shown for elastic proton scattering from 16 at 500 MeV 
laboratory energy. The solid line represent the calculation performed with a first-order full-folding 
optical potential as described in Section II. The dashed line represents a calculation where the 
M0ller factor is evaluated for the fixed angle = 90°, whereas for the dotted line the M0ller factor 
was omitted altogether. All calculations are based on the DH density and the Bonn model D52 
t-matrix. 

FIG. 7. The angular distribution of the differential cross-section (^), analyzing power (A y ) 
and spin rotation function (Q) are shown for elastic proton scattering from 40 Ca at 100 MeV 
laboratory energy. The solid line represent the calculation performed with a first-order full-folding 
optical potential based on the DH density [16] and the full Bonn model f22]| , the dashed curve is 



based on the factorized, off-shell approximations. The data are taken from Ref. [31] 



FIG. 8. Same as Fig. 7, except that the HFB model [17] is employed for the density as well 
as the mean field force. 
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FIG. 9. The angular distribution of the differential cross-section analyzing power (A y ) and 
spin rotation function (Q) are shown for elastic proton scattering from 90 Zr at 80 MeV laboratory 
energy. The solid line represent the calculation performed with a first-order full-folding optical 



potential based on the DH density |16| and the full Bonn model [22|, the dashed curve is based on 



the factorized, off-shell ""tp 1 approximations. The dotted line represents a full-folding calculations 



based on the HFB model [17]. The data are taken from Ref. [32]. 



FIG. 10. The angular distribution of the differential cross-section (^), analyzing power (^4^) 
and spin rotation function (Q) are shown for elastic proton scattering from 208 Pb at 65 MeV 
laboratory energy. The solid line represent the calculation performed with a first-order full-folding 
optical potential based on the DH density [16] and the full Bonn model |2^| , the dashed curve is 



based on the factorized, off-shell 'ip' approximations. The data are taken from Ref. [33]. 



FIG. 11. Real and imaginary part of the on-shell value of the optical potential as function 
of the orbital angular momentum L for scattering from 40 Ca at 200 MeV laboratory energy. (+) 
denotes the potential for J = L + |, whereas (-) stand for J = L — i. The full- folding (solid) and 
factorized, off-shell (dashed) calculations are based on the full Bonn model and the DH density. 

FIG. 12. Real and imaginary part of the on-shell value of the proton-nucleus t-matrix (coulomb 
contributions omitted) as function of the orbital angular momentum L for scattering from 40 Ca 
at 200 MeV laboratory energy. (+) denotes the potential for J = L + ^, whereas (-) stand for 
J = L — i. The full-folding (solid) and factorized, off-shell (dashed) calculations are based on the 
full Bonn model and the DH density. 

FIG. 13. The angular distribution of the differential cross-section (^), analyzing power (A y ) 
and spin rotation function (Q) are shown for elastic proton scattering from 40 Ca at 200 MeV 
laboratory energy. The solid line represent the calculation performed with a first-order full-folding 



optical potential based on the DH density [16] and the full Bonn model 122], the dashed curve is 



based on the factorized, off-shell 'ip' approximation. The data are taken from Ref. [34 

FIG. 14. Same as Fig. 11, except for 208 Pb. 
FIG. 15. Same as Fig. 12, except for 208 Pb. 



FIG. 16. Same as Fig. 13, except for Pb. The data are taken from Ref. [35] 
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FIG. 17. The neutron-nucleus total cross-sections for scattering from 12 C, 16 0, 28 Si, 40 Ca, 90 Zr, 
and 208 Pb are shown as a function of the incident neutron kinetic energy. The dotted line represents 
the data taken from Ref. [28,2{|. The solid diamonds correspond to the full- folding calculations 
using the full Bonn NN t-matrix [22| and the DH model fl6| for the density. The open circles 
correspond to the factorized, off-shell 'i// approximation. The calculations for energies smaller and 
equal 200 MeV include the propagator modification due to the DH mean field. 
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